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We present a quantum theory of low-lying excitations in a trapped Bose condensate with finite 
particle numbers. We find that even at zero temperature condensate number fluctuations and/or 
■ fluctuations of the excitation frequency due to quantum uncertainties of the mode occupation lead 

^\ ' to a collapse of the collective modes due to dephasing. Coherent revivals of the collective excitations 

, are predicted on a much longer timescale. Depletion of collective modes due to second-harmonic 

t-H ■ generation is discussed. 
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I. INTRODUCTION 

> ■ 

In order to investigate the properties of the recently created Bose condensates of alkali atoms in magnetic traps 
[Q-Q the excitation of low-lying collective modes by periodic variations of the trap potential has turned out to be a 
very effective tool. The successful implementation of this idea has produced results for the collective mode spectrum 
of the trapped condensate [j|,p [which are in very good agreement with theoretical predictions based on solutions of 
the Gross-Pitaevski equation [7|-Q i.e. on classical mean field theory at zero temperature. In the experiments, in 
addition to the mode frequencies, the decay of the collective excitations in real time could be observed. At present 
there exists no theory which describes the damping of excitations at finite temperature. In this paper we focus solely 
on the decay of the collective exciations at very low (effectively zero) temperature. In the very low temperature 
regime it is difficult to imagine any truly dissipative mechanism for the observed decay, as long as the trapped system 
I ■ can be considered closed. The usual dissipative mechanisms in a homogeneous Bose condensate depend on the 
' O \ continuous mode spectrum found there. 

It has recently been shown |Q that nondissipative interactions give rise to collapses and revivals of the macroscopic 
wavefunction for small atomic condensates, analogous to those predicted |ll| and observed |14| for a single mode field 

• • . and a two-level atom. Here we extend this idea to show that very similar mechanisms may lead to the collapse 
(apparent damping) of collective excitations in a finite Bose condensate with a discrete mode spectrum. For times 
much longer than the observation time in present experiments, this collapse is predicted to be reversed in 'revivals' 

^ ■ of the collective excitations. 

We shall investigate two possible mechanisms for the collapse. The first mechanism, discussed in Section II, is 
based on atom-number uncertainty in the condensate. It therefore has the same physical origin as the collapse of the 
macroscopic wave function However, it will turn out that this mechanism is much less effective for collective 

modes. The second mechanism, described in Section III, is based on a potential dependence of the mode frequency on 
the mode occupation. Given such a dependence, quantum uncertainties in the occupation of a mode must lead to its 
collapse. A two-mode model with this property has been discussed recently by Kuklov et. al. jl5), while Pitaevskii Jl]]] 
also recently discussed this mechanism in general for nonlinear oscillators and applied it to a special collective mode in 
trapped condensates. Here we wish to present a theory of the nonlinear self-coupling of the collective hydrodynamic 
modes in trapped Bose-condensates and examine the conditions under which the nonlinear coupling can become 
effective. 

We emphasize that while we propose these collapse-mechanisms as possible explanations for the damping of exci- 
tations at zero or very low temperatures, in current experiments there appears to be a substantial finite temperature 
contribution which is likely to produce the dominant contribution to the currently observed damping times [18). 
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II. COLLAPSE DUE TO ATOM-NUMBER UNCERTAINTY IN THE CONDENSATE 



The mechanisms we present are intrinsically quantum mechanical and, in principle, occur in any finite Bose con- 
densed system, in which the spontaneous symmetry breaking associated with the infinite system cannot occur [ jl9| . 
For this reason the usual Bogoliubov analysis of the low-lying excitations of an infinitely extended Bose gas, which is 
based on the assumption of a spontaneously broken gauge symmetry, is not directly applicable. It must be modified 
by eliminating the underlying assumption of broken gauge symmetry before the collapse mechanism proposed here 
can be described consistently. Here we first present this modification of the Bogoliubov theory, then examine the 
collapse and revival of the collective excitations. 

The grand canonical Hamiltonian of the system is [|20fl 



H = d 



v3„ 



+-2 T J 

— Vft-Vi> + (V(r) + SV(r, t) - u) $U + -?-^ 2 ^ 2 
2m 2 



(1) 



where m is the atomic mass, fi the chemical potential, V(r) is the static trap potential, SV(r, t) its modulation, and 
Uq = Air?i 2 a/m is proportional to the s-wave scattering length a. The Bogoliubov aproach to the low-lying excitations 
is based on decomposing the boson field annihilation operator -0(r) (and its adjoint) as f2l"| , [22f 

V3(r) = (V3(r))+^(r), (2) 

and then approximating ([j]) by a quadratic form in the operators Sip(r) and dffi(r). This approach is successful if the 
expectation value (V>(r)) on the right-hand side of (||) is nonzero, which is the case in an infinitely extended system if 
the U(l) gauge symmetry is spontaneously broken, or in a finite system if the phase of the condensate is established 
with respect to a reference, e.g. by a measurement of its phase relative to an infinitely-extended Bose-condensed 
system filifl , which itself has a well-defined phase. However, in the experiments on collective excitations in trapped 
Bose condensates referred to above a phase preparation of the condensate is not made, so that (V>(r)) = on the 
right-hand side of (^) and this decomposition becomes useless. In fact, the same kind of problem appears for a laser 
far above threshold, where Bose condensation of photons in the laser mode occurs, but the phase is not fixed and 
undergoes a diffusion process p3| . As is done in that case, it is necessary here to replace (||) by [] 

&r) = e^^p a (r)+6p(r)+Mr) . (3) 

The decomposition (3) was used by Popov p3] for spatially homogeneous Bose condensed systems. Here Po( r ) is a c_ 
number and Sp(r) and (j>(r) are operators with commutation relations [Sp(r), e l *( r '] = — e^^S^ (r — r') and together 
they provide a quantum mechanical description of the condensate and its collective excitations at long wavelengths 
pjj . It is convenient to decompose 6p(r) further into 

5p{r) = 8a(r) + 5n/Q , (4) 

where Sn — J d 3 rSp(r, t) is canonically conjugate to the spatially averaged phase variable 4> = / d 3 r<fi(r) according 

to the commutaion relation [Sn, <fi] = i and therefore commutes with V</> and has discrete, not necessarily integer 
eigenvalues with unit spacing. It describes number fluctuations in the condensate. In the following it is useful to 
introduce the abbreviation p(r) = po(r) + 5h/£l so that po + Sp = p + Sa. Note that p(r) also commutes with V(f>. The 
operator "0i(r) appearing in Eq.(|^) describes short- wavelength components of $(r) whose elimination renormalizes 
the coefficients of the effective long- wavelength theory for 4>(r), 5p(r) |24f| . and which provide a particle reservoir for 
the condensate even in the case T = to which we confine ourselves here. In the following we shall assume that 
the elimination of ipi has already been performed. The ansatz (^) is used in the Hamiltonian (|l|) to determine the 
Heisenberg equation of motion for Sp and <fi. We obtain formally (i.e. without paying attention to the fact that 
operator products at equal points in space are only well defined after chosing an apropriate regularization procedure; 
for our present purposes this is sufficient, since we are only interested in the quantized fields Sp, 6(f) either in their 



1 The decomposition of t/j and ^ into amplitude and phase is beset by well-known problems, which arise due to the existence 
of the vacuum state satisfying ij)\Q) = 0. However, since No = J d 3 rpo(r) S> 1 the probability amplitude for the vacuum state 
1 0) is negligibly small. 
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free-field regime, like in the present section, or in the case where all modes except one are in their vacuum states, as 
in the following section), 



5p{v, t) = -(fi/m) V- Mr) + Sp(r, t)) 1 ' 2 V0(r, t) (p (r) + Sp(v, <)) 



,1/2 



(5) 



$(r, t) = -~(U (p (r) + Sp(r, t)) + V + SV(r, t) - p) 

~[(fio(r) + Sp(r, t)y 1/2 (V0) 2 (po(r) + Sp(r, *)) 
4m 



It follows from 



that 



1/2 



+ ( Po (r) + Sp(r, t))- 1 ' 2 (V^) 2 ( Po (r) + Sp(r, t))- 1 ' 2 ] 

(po(r) + ,5/3(r ' t)rl/2 (v)2 (/9o(r) + ,5/3(r ' <))1/2 



(6) 



5n 



d s rSp = . 



(7) 



i.e. there is no restoring force on Sh, which therefore need not necessarily be small. Therefore we combine this 
quantity with po to extract from (o) the equation which holds to zeroth-order in the small quantities Vcf> and Scr 



pp(r) 



4m 



(Vp(r) - ip(r)- 1 (Vp(r)) 2 \ + Vp(v) + U p(rf 



(8) 



This is one component of the Gross-Pitaevski equation and identifies p(r) as the operator of the density of the 
condensate, as well as defining the operator of the chemical potential p — p(N + 8h) as a function of the number 
operator N + Sn — J d 3 rp(r) in the condensate. Next we linearise with respect to the small quantities 5a(r,t) and 
V(/>(r) to obtain 



5&(r,t) = V- 

m 



p(r)V(f>(v,t) 



(9) 



±U 5*(r,t) + ^P(r)" 1/2 



(v 2 - p{v)- l/2 {V 2 p{r) 1/2 )} p(r)-^ 2 5a(r, t) - ^SV(r, t) . (10) 



Eliminating V0(r) from (^) and (JT^) we then get the wave equation for the low-lying excitations, 

S'a{r, t) = -V- [p(r)V (U Sa(v, t) + SV(r, t))} 
m 

p(v)V \p(r)- 1/2 fv 2 - piT)- 1 ' 2 ^ 2 ^) 1 ' 2 )) p(r)-^6a(r,t) 



Am 2 



(11) 



Neglecting the gradient term of fourth-order compared to those of second-order, and also neglecting Sn compared 
to po and replacing 5a by the classical density fluctuation, we recover Stringari's wave equation jTJ] for the collective 
excitations of a trapped Bose condensate. We shall in the following also neglect the higher-order gradient terms, 
but shall keep the number fluctuation operator Sn, which is implicit in p, and examine its consequences in Eq. (|T]]). 
Since the short-wavelength non-condensate components of the system act as a reservoir even at T = these number 
fluctuations may reasonably be taken as Poissonian. We can therefore assume that the initial state of the system 
is given by a pure state or a mixture with a Gaussian distribution over eigenstates of Sh, with zero mean particle 
fluctuation and mean square deviation given by ((AiVo) 2 ) = No, the average number of particles in the condensate. 
This corresponds to a Poissonian distribution of the total particle number in the condensate Sh + iVo around its 
mean value Nq ^> 1. We emphasize that due to the entanglement of the condensate with the other states of the 
system, which have been eliminated, the use of a mixture of Sh eigenstates applies even to a single realization of the 
experiment, and not only to an ensemble of realizations corresponding to identically prepared experiments. For each 
eigenvalue Sn of Sh we have in principle to determine the solution of Eq. (||) and of Eq. ( p"l| ) corresponding to a given 
mode. This gives the normal modes of the density oscillations, after SV(r, t) has been switched off, and determines 
their frequencies ui v as a function ujq^{Nq + Sh), where u>o v {Nq) is the collective mode frequency determined by the 
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mean- field theory fi^-|lO|]. In the quantum ensemble defined by the initial state only eigenvalues Sn of Sn which are 
very small compared to Nq are contained with appreciable weight. Therefore it is sufficient to expand w„ to first-order 
in Sn 

lu u = W0i/(1 + IvSn) , (12) 

with 

a quantity determined by the solutions of the mean-field theory [fz|— p~p|| . 

We can now discuss the collapse and revivals of the collective excitations by evaluating the quantum ensemble 
average of the density oscillation (6a u (t)) for a given mode, assuming that it is coherently excited, e.g. by modulating 
the trap at the required frequency at times prior to t = 0, while for t > the mode is left to evolve freely. We obtain 
at times < t < (-y^wou)^ 

(Sa„(t)) = exp(-ijVo7^t 2 ) f (6& v (0)) cos(u 0v t) + s i n (oj 0v t) ) , (14) 



It can be checked that the corrections of order (Sn) 2 to (12) make a negligible contribution in (|Tj). Equation (|l 



shows that the excitation decays by dephasing on a timescale 



r c = (y/No/2h»\'*»>)~ 1 , (15) 

in the form of a Gaussian collapse. This collapse occurs even at a temperature T = 0, and even if the condensate 
never has a well-defined overall phase. This collapse is therefore different from (and turns out to be less effective 
than) the collapse of the macroscopic wave function due to number-fluctuations [|12| , which occurs if a well defined 
phase of the condensate (with respect to some phase standard, like another condensate) is prepared e.g. at t = 0. On 
timescales t > (jvUiou) -1 the discreteness of the spectrum of Sn manifests itself and we obtain revivals of (8p v (t)) at 
the times t — mr y lr y u uJov . 

To estimate the order of magnitude of the effect and to examine its accessibility to observation we have to obtain 
an estimate of 7„ e.g. for the observed m = modes in the experimentally realized condensates. We can achieve this 
goal by using some analytical results due to Stringari H . It was shown in Joj that for N — > oo the mode frequencies 
become independent of No. It was also shown there by using sum-rule arguments that for finite No the loov for the 
low-lying states can be represented in the form u>qv(No) = wo(°°)(l + c u Eki n / Eho), where the c„ are parameters of 
order 1 which depend on the trap geometry and the particular mode, and Eki n /Eh is the ratio of the kinetic energy 
and the harmonic trap energy in the ground state. The No dependence of this ratio in the limit of large Nq of interest 
here can be estimated by using the Thomas-Fermi approximation |25f| as E^i n j E) lo = b v {No a l \/^/niuJo)~ 4 ^ 5 with ujo 
the harmonic trap frequency and a proportionality factor b v of order 1 which is asymptotically independent of No but 
dependent on the trap geometry. For the coefficient 7„ we obtain finally 

lu = -%c„N- 9/5 (a/d Q r 4 / 5 , (16) 
5 



with do = \/S/mwo, giving a collapse time 

/V 13/10 

t c = (5V2/4b^)(a/d )^ 5 ^ , (17) 

which is larger than the collapse time for the macroscopic wave function obtained in [Q roughly by a factor of fx/u>o v - 
For the two experiments in which decay times of collective excitations have been measured, the following parameters 
can be estimated: In the experiment of [^), the lifetime of the m = mode (with a frequency of vq v — ^o^/2vr = 1.84tv 
where the radial trap frequency has the value v r = 132Hz) for a condensate containing No = 4500 rubidium-87 atoms 
was measured to be 110 ± 25ms. For this case the following numbers apply: ujovI^ = l87Hz,No = 4500, a = 
52A, U>o = LJr (cj z /w r ) 1 / 3 ,a/ v /^Wo = 10~ 2 . This gives a collapse time of r c = 850ms, which is longer than the 
observed decay time, but this is to be expected, as the latter was measured in a regime where finite temperature effect 
are not negligible |Tq |. 

In the experiment of the lifetime of the 30Hz collective excitation of a condensate of iVo = 5 x 10 6 sodium atoms 
was measured to be 250ms. Due to the much larger number of atoms used in this experiment the mode frequencies 
loov become insensitive to the dispersion in the number of particles in the condensate || , and the collapse due to the 
mechanism under discussion here occurs only on time-scales very long compared to the observed damping time. 
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III. COLLAPSE DUE TO NONLINEARITY OF THE COLLECTIVE MODES 



In this section we neglect any uncertainty in the number of particles in the condensate. However, we no longer 
linearize the Heisenberg equations of motion for 5p and 8<j>. Instead we consider the case where only a single mode of 
the collective excitations is appreciably excited and ask for the nonlinear corrections to the dynamics of this mode. 
This nonlinearity will also give rise to a dephasing and collapse of the collective mode amplitude. This is in the same 
spirit as the recent work of Kuklov et. al. and Pitaevskii ]l7jj, but differs from the former in that we don't treat 
a model but rather start from the full microscopic approach, and is more general than the latter because the analysis 
is not confined to a special symmetric collective mode. 

The analysis starts with the nonlinear Heisenberg equations of motion Eqs. (||) and ([}]) but with the density 
gradient-terms neglected. These can be derived from the Hamiltonian 

(18) 



Questions of operator ordering are not important at this stage. In the end we will choose normal ordering in the 
mode operators; this is the meaning of the notation {. . We introduce a mode expansion using the modes of the 
linearised hydrodynamics, 



w = i ?'( 



/,J V„(r) tt „(f)-h.c. , (19) 



2U, 



{ ° V-F„(r)a v (t)+h.c. J , (20) 



2hu 



where ^ is the mode sum not including the spatially and temporally constant part; the constant part is already 
included in p = po +Sn/Q. Here u u is the frequency associated with the linear mode function F v (r), so that the linear 
hydrodynamics H implies a v {t) — a v eT lUvt . 
The Hamiltonian now becomes 

ff = £W^ + ftWic« A) c*^^ , (21) 

V VK.X ' 

with C*- 1 ) and symmetric in the parenthesised indices and given by 

,(i) _ ifi /{To"/ (i) (!) (i) 



with 



and 



where 



C) x > ^ = —\l—( D (1 } + D (1 ) x . + ,) , (22) 



D SU = I d " xFv (VFk • VFx) ' (23) 



(2) _ ^_ JUo ( 2D (2) (3) D (3) D (3) „(3) \ (24] 



v ' y uj k lu\ j y ui k oj\ j 

The resulting equations of motion in the interaction representation are 
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We now assume that only a single mode [i is externally excited; all other modes are only excited via their coup ling 
to the mode \i. This implies that can be considered "large" compared to all other mode operators. Using Eq. (|2q ) 
in this approximation, we have for v = /i 

a, * -iJ2> {2C^ K /^^ a y K + cW K) e-^a»a K + C^^afa + C e<(H ""M««} - ( 27 ) 
and for v = k ^ fi 

- ~i {og^^W + 5^)^"-^ J W + Cg^e^aj,^} . (28) 

In general, the oscillatory terms in this will be changing much more rapidly than a M in the interaction picture. The 
only exception to this is when we have a second harmonic resonance, lu k = 20^; this case needs separate treatment, 
and is discussed below. There is no contribution from the apparent resonance when oj k = 0, since this mode is 
explicitly excluded from the sum; in fact, the coefficient C^)q would be zero anyway. Thus we may solve Eq. fl2S| ) 
for a K treating a M as approximately constant, to find 

c (l)* c (2) c (2)* } 

j K + 2w M M 2uj k - 2lj m w k m 



We now subsitute this back into Eq. (g7|) for a M , and keep only the nonoscillatory terms, to obtain 

A/ — i kcrv^ i 



(30) 



where 



.(1) 1 2 |^( 2 ) 1 2 , l/-( 2 ) " 



The form of the matrix elements (22)-(24) implies selection rules for the modes k contributing to the sum in Eq.(31). 
For axially symmetric trap potentials with inversion symmetry only modes contribute with positive parity and az- 
imuthal quantum numbers m K satisfying m K = 2ra^ or m K — (via C ? % ^ 0) or m K = -2m,, ( via C9^s ^ 0). 

The effective self-coupling coefficient k experimentally manifests itself by an energy dependence of the observed 
mode frequency according to uj^E^) = w p — (n/huo^E^. Having derived (pi|), the collapse and revival of the mode 
/i follow immediately: If at t = the mode is excited e.g. in a coherent state with amplitude A then for t > its 
amplitude changes according to 

(A|a^|A) = Aexp [-|A| 2 (1 - cos(/tf))] [cos(|A| 2 sin(>rf)) - isin(|A| 2 sin(/ct))] . (32) 

which for times \n\t « 1 collapses with the new collapse time r c = (|Ak|) _1 according to (A|a^|A) w Aexp(— \{t/t c ) 2 
but is revived at revival times t r = 2nirK~ 1 for integers n > 1. 

We can estimate the order of magnitude of the effective coupling constant k, in particular its scaling with the 
system parameters and the number of atoms in the condensate, by using the Thomas-Fermi approximation 

K J^^L(±.Y^ oN -^(^) V \ (33) 
uj m^w uj q \r TF J \ a J 

which gives for the collapse time the estimate 

2/5 



In this result, derived for an initially excited coherent state, we may replace the coherent amplitude A by the variance 
of the excited quantum number to generalize it for an arbitrary initially excited state. The scaling of k with the 
system parameters is consistent with Pitaevskii Jl7t who derived the coupling coefficient for a special mode. The 
scaling of the coefficient with an inverse power of Nq tends to make it small in most cases. Putting in the numbers 
for the two experiments, assuming the excited quantum number N e to be 1% of the total particle number we obtain 



G 



Rb: N e = 45, A = ^45, t c = 2.2s 
Na:N e = 5- 10 4 , A = 10 2 Vb, t c = 1.47 • 10 3 s. 



(35) 



which is longer than the observed decay 

Relatively larger values are obtained for k (and shorter ones for r c ) if there is a collective mode near the second 
harmonic of the externally excited mode. Second-harmonic generation has recently been seen in numerical simulations 
of the Gross-Pitaevskii equation J27J . 

We shall briefly discuss this point for the case of an anisotropic axially symmetric trap whose collective mode 
frequencies and eigenfunctions are known as a function of the anisotropy parameter [3 — lo z /loq in the Thomas- Fermi 
and hydrodynamic limit ( [p8[), where w 2 ,wo are the axial and radial trap frequencies, respectively. The modes 
are in this case labelled by three quantum numbers (n, j, m) (see p8|]). We are interested in the case where for two 
different sets of their values we have second harmonic resonance such that io v = 2c^, and where the relevant matrix 

clcn lent Cr/ * with v = [n,j,m),(i = (n,j,m) does not vanish. The latter condition implies certain selection rules, 
which for the case at hand read: i) The mode function F u (r) of the second harmonic mode in cylindrical coordinates 

(2) (3) 

p, z, <p must be even in z, and ii) at least one of the two conditions fn — 2m (for D^/ ■. ^ 0, ^ 0) or fn = 0(for 
D(3) ^ 0) must be satisfied. 

As an example we consider the mode p = (n = 2,j = 0, m = 2) which has uj^ = \/2ujo, -P^( r ) — N^p 2 exp(2i<p) 
where iV M is a normalization factor, and we shall fix the phases of the mode functions by chosing these nor- 
malization factors always positive. Then the two modes vl = (n = 2,j = l,m = 0) with lo^ x = (3/3 2 /2 + 
2 + y/(2 - 3/3 2 /2) 2 + 2( 3 2 )w%, F„i(r) = Nui(-Ap 2 /3 - 16z 2 /3 + 1) and u 2 = (n = 2,j = 0,m = 4) with 
uj 2 u2 = (3/3 2 /2 + 10 - ^(6- 3/3 2 /2) 2 + 10/3 2 )^, F v2 (r) = N v2 (-lQz 2 + 2p 2 - l) /9 4 exp(4^) both satisfy the selec- 
tion rules. For the special value /3 2 = (ui z /uiq) 2 = 16/7 it turns out that both modes vl,v2 are degenerate and 
in second harmonic resonance with the mode at frequency lu^ = v2woi i-e. Wvi = <j-V2 = 2\/2ujq. (A resonance 
for this value of (3 2 has recently also been noted in |2^]). The relevant matrix elements in this case turn out to be 
C i%n) = 0.562A ( 7 7/10 (d /a) 1 /5 Wo and = 0.332A~ 7/10 (do/a) 1/5 wo For such a case of exact second harmonic 

resonance, the previous treatment breaks down. The differential equations for the two, or, in the present case even 
three, resonantly coupled modes must instead be considered together. Neglecting couplings to all other modes, we 
have for w vi — 2lu^ 



*•* = ~l c i 2 L)< > ( 37 ) 

which are the coupled mode equations for second harmonic or subharmonic generation. It can be seen that only the 
linear combination a„ — const. 2_ li C^J> *.a V i of second harmonic modes couples to the fundamental mode p, while the 

linear combination a„ = const\C^\sa. v \ — C^^^a^] is not generated. Choosing the constants (const = | _1 ) 
to normalize these linear combinations, Eqs.(36),(37) reduce to 



= -i||cSLl«t«- (38) 



-IT (2) In- 2 



(39) 



where we defined \C 



(2) 



\C. 



(2) 

«2(/i/i) 



j 2 ] 1 / 2 , which takes the value |C^' 



(2) 
(tili) I 



0.6537V, 



-7/10 







(d /a 



.1/5. 



loq m 



our example. These equations have well-known exact solutions 1 30 1 in terms of Jacobi elliptic functions, that oscillate 
with a frequency 



An order-of-magnitude estimate of this is 



m 2 toor^ F 



— » r - 7/10 , -v. 



1/5 



(40) 



(41) 
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which is much larger than the nonresonant rate K^rT^. The initial transfer of energy from the fundamental mode to 
the second harmonic may be viewed as a collapse with a collapse time r c ~ fio- 

The resonance phenomenon found here might help to explain why in the experiments some modes show amplitude 
dependence of the frequency, while others do not. 

IV. CONCLUSIONS 

We have investigated two different quantum mechanisms which may lead to a dephasing of collective modes in 
trapped Bose-Einstein condensates even at zero temperature. Both rely on the possibility that the system may be in 
a linear superposition of states each of which has a slightly different frequency for the collective mode. In the first 
case investigated here the linear superposition combines states with different numbers of atoms in the condensate. 
This may occur since the non-condensate atoms, which are always present due to many-body interactions pi] , act 
as a particle reservoir for the condensate. This case was earlier shown to give rise to collapses and revivals of the 
macroscopic wave-function once the latter is prepared at a given time with a well-defined (relative) phase, e.g. 
by a measurement. Here we found that a dephasing effect occurs also for the collective mode but on a different time 
scale which is typically much longer for large condensates. In principle, the collapse of the macroscopic wave-function, 
if it was initially prepared with a well-defined phase, and of the collective mode amplitude may occur independently 
and simultaneously, but their simultaneous observation would of course be much more difficult than the observation 
of each collapse individually. While the dimensionless collapse time of the macroscopic wave function for large Nq 
scales like 

uqt c r~j (do I a) 2 / 5 N^ 5 / (ANq) (macroscopic wave function) (42) 
the collapse time of the collective mode based on this mechanism scales as 

uqt c ~ (do / a)~ 4 / 5 N^ 5 / (ANq) (collective mode due to atom-number uncertainty). (43) 

1 /2 

which corresponds to Eq.(17), which was written for (AN ) = N > . 

In the second case investigated by us the linear superposition is one of different quantum number states of the 
excited collective mode, which will typically be in a coherent state immediately after its excitation from its vacuum 
state. We have calculated the effective nonlinearity of the excited mode due to its coupling to the other modes. This 
nonlinearity gives rise to a dispersion of the collective mode frequency within the linear superposition of number 
states. The collapse time obtained from this mechanism we find (in Eq.(34)) to scale like 

wot c ~ (do/a)^ 2 ^ 5 N^ 5 /(An) (collective mode due to nonlinearity) (44) 

For an initially excited coherent state the variance of the quantum number is related to the average mode-energy 
by (An) = y/ E^/Tiuj^. This scaling is consistent with a result obtained by Pitaevskii ]l7|| . We have also found that 
the nonlinearity may be strongly enhanced if there is a mode at or close to a resonance with the second harmonic of 
the excited mode. In the case of such a resonance the transfer of energy from the fundamental mode to the second 
harmonic for short times looks similar to a collapse of the fundamental mode, which, according to Eq.(41), occurs on 
a time-scale scaling like 

uqt c ~ (do/a) -1 / 5 A^ 10 /^/ E^/hu!^ (collective mode due to second-harmonic). (45) 

Due to the scaling of the collapse times with positive powers of the atom number the collapse by one of the mechanisms 
we have examined could best be observed in comparatively small condensates. At presently achieved temperatures 
the observed damping of the collective mode occurs at a higher rate than the collapse rates obtained here. However 
the observed damping rates were found to be strongly temperature dependent, decreasing rapidly with decreasing 
temperatures jj^]. Therefore the collapse we have investigated here is still masked in present experiments by finite 
temperature effects, but is predicted to reveal itself when experiments are pushed to smaller temperatures. 
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